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ABSTRACT 


To  make  a  drug-testing  program  suc¬ 
cessful  and  to  minimize  the  cost  of  the 
program,  the  minimum  number  of  tests  that 
must  be  given  in  a  specified  period  to  identify 
a  fixed  percentage  of  drug  users  must  be 
determined.  This  memorandum  presents  a 
Markov  model  that  can  be  used  to  determine 
the  number  of  tests  that  should  be  given.  In 
addition,  three  applications  of  the  model, 
showing  how  it  can  be  used  to  analyze  the 
drug-user  population,  are  presented. 
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INTRODUCTION 


The  ability  of  urine  tests  currently  used  to  identify  drug  users  depends 
on  the  level  of  the  drug  present  in  the  subject’s  body  when  the  test  is 
administered.  The  level  of  drug  in  the  body  decreases  over  time  after  the 
substance  is  used.  Thus,  how  often  an  individual  uses  a  drug  and  the  time 
lapse  between  drug  use  and  a  drug  test  are  important  variables  in  predicting 
how  successful  a  testing  program  will  be  in  identifying  users.  To  minimize  the 
cost  of  a  program  of  drug  testing,  the  minimum  number  of  tests  that  must  be 
given  over  a  specified  period  in  order  to  identify  a  set  percentage  of  users  in  a 
population  must  be  determined. 

This  memorandum  introduces  a  Markov  model  developed  to  determine 
the  minimtim  number  of  tests  that  should  be  given.  Some  of  the  assumptions 
made  in  developing  the  model  impose  certain  limitations  on  its  use.  In 
particular,  the  model  assumes  that  a  drug  test  given  by  the  Navy  will  have  no 
"false  positives”;  that  is,  the  test  will  not  indicate  drug  use  if  given  to  a  drug- 
free  individual. 

Applications  of  the  model  in  analyzing  the  drug-user  population  are  also 
presented.  The  model  is  a  generic  one,  and  the  examples  of  its  applications 
presented  here  are  meant  to  illustrate  its  potential  uses.  Some  of  the 
assumptions  made  in  the  examples  about  the  distribution  of  drug  users  in 
certain  populations  may  not  conform  to  reality.  Most  notable  of  these  assump¬ 
tions  is  that  distribution  of  drug  users  among  new  recruits  is  identical  to  the 
distribution  of  drug  users  in  the  Navy  at  any  given  time.  In  any  actual 
application,  the  model  should  be  modified  to  reflect,  to  the  extent  possible, 
characteristics  of  the  drug  tests  used  and  the  populations  being  tested. 

BASIC  ASSUMPTIONS  AND  NOTATIONS 

The  Navy’s  current  drug  tests  can  detect  a  user  who  has  taken  the  drug 
no  more  than  a  specified  number  of  days  before  the  test  is  administered.  For 
purposes  of  this  discussion,  it  is  assumed  that  the  test  used  will  detect  a  user  if 
he  has  used  the  drug  in  the  last  n  days.  Also,  the  probability  of  detection  is 
assumed  to  be  a  function  of  the  number  of  days  since  the  drug  was  last  used. 
In  particular,  d(i),  i  =  1, ...,  n,  will  denote  the  probability  that  the  drug  will  be 
detected  by  the  test  if  the  person  being  tested  last  used  the  drug  i  days  ago. 


The  distribution  of  drug  users  within  the  group  to  be  tested  is  also 
assumed  to  be  known.  In  particular,  it  is  assumed  that  users  have  been 
classified  according  to  their  usage  rates— for  example,  once  a  day,  once  a 
week,  or  once  a  year.  The  number  of  distinct  classes  of  users  is  assiuned  to  be 
m  +  1,  denoted  by  0,  1,  2,  ...,  m,  and  class  0  corresponds  to  the  class  of 
nonusers.  The  probability  that  an  individual  in  class jf,  j  =  0, ...,  m,  will  use 
the  drug  on  an  arbitrarily  selected  day  will  be  denoted  by  p(j).  The  probability 
that  a  person  selected  at  random  will  be  in  classj  will  be  denoted  by  rij),  j  =  0, 
1, ...,  m.  Note  that  the  probability  of  a  person  being  in  class  =  1,  2, ...,  m, 
given  that  the  person  is  a  drug  user,  is 


r(l)+  r  (2)  +  ...  +  r  (m) 


THE  MARKOV  MODEL 

A  finite-state  Markov  chain  is  used  to  estimate  the  probability  of 
detection  when  the  drug  test  is  given  at  random  to  a  person  in  one  of  the  usage 
classes.  The  states  of  the  Markov  chain  correspond  to  the  number  of  days 
since  the  drug  was  last  used.  In  particular,  a  person  in  class  j,  j  -  0, ...,  m,  can 
be  in  n  +  1  states,  1, 2, n  +  1,  where  state  i,  i  =  1, ...,  n,  corresponds  to  use 
of  the  drug  exactly  i  days  ago,  and  state  n  +  1  corresponds  to  last  use  of  the 
drug  at  least  n  +  1  days  ago  or  never  having  used  the  drug. 

If  an  individual  from  class  j  is  in  state  i,  i  =  1,  ...,  n,  today,  the 
probability  that  he  will  be  in  state  1  tomorrow  is  p(jh  The  probability  that  he 
will  be  in  state  i  +  1  tomorrow  is  1  -  p(j).  The  probability  of  his  entering  any 
state  other  than  1  or  i  +  1  is  0.  If  he  is  in  state  n  +  1,  the  probability  of 
entering  state  1  tomorrow  is  pO),  and  the  probability  of  being  in  state  n  +  1 
tomorrow  is  1  -  p(j).  The  probability  of  going  to  a  state  other  than  1  or  n  1 
is  0.  If  pO)  is  not  zero,  the  following  Markov  matrix  results. 


Beginning 

Ending  state 

state 

1 

2 

3  n- 1 

n 

n  +  1 

1 

P(j) 

1  -p(j) 

0  0 

0 

0 

2 

P(j) 

0 

1-P(j)  0 

0 

0 

■ 

I 


n-  1 

p(i) 

0 

0 

0 

>-p(i) 

0 

n 

p(j) 

0 

0 

0 

0 

1  -p(j) 

n  ♦  1 

p(j) 

0 

0 

0 

0 

1  -p{j) 

This  Markov  matrix  has  the  following  steady-state  probabilities 
associated  with  it: 


x(i)  =  pO) 
x(2)  =  p(/)ll  -p(/)) 

x(3)  =  pO‘)11-p(7)1^ 
x(4)  =  p(/')tl-p(/)r’ 


x(n-l)=pO-Xl-pO‘)l‘'’"^’ 
x(n)  =  p011-p(/)r"^’ 
x(n+l)  =  [l-pO')j" 


Uj  =  0,  then  p(j)  =  0,  x{i)  =  0  for  i  =  1, n,  and  x(n  -I-  1)  =  1.  Hence, 
when  a  test  is  administered  to  an  individual  from  class  ;  =  1,  m,  the 


probability  of  detecting  drugs  in  this  individual  is 

D(j)  =  d(l)x(l)  +  ...  +  d(n— 1)  +  dBin)x{n) 

D(0)  =  0  . 

The  probability  that  an  arbitrarily  chosen  user  is  detected  is  given  by 

Du  =  s(l)i)(l)  +  s(2)D(2)  +  ...  +  s(m— l)D(m  — 1)  +  s(m)D(m)  . 

Du  can  be  interpreted  as  the  detection  rate  among  users.  The  probability  that 
a  person  selected  at  random  will  show  drug  usage  on  a  test  (i.e.,  the  rate  of 
detection  over  the  entire  population)  is  given  by 

D  =  r(l)D(l)  +  K2)D(2)  +  ...  +  dm-l)D(m-l)+  tim)D{m)  . 

MINIMAL-COST  TESTING 

The  model  developed  above  can  be  used  to  determine  the  smallest 
number  of  tests  that  can  be  given  to  each  individual  in  the  population  to 
guarantee  that  a  certain  percentage  of  drug  users  will  be  identified.  Let  t. 
0  <t  <  1,  be  the  decimal  equivalent  of  the  percentage  of  actual  drug  users  we 
wish  to  detect.  Assuming  that  the  tests  will  be  administered  to  individuals 
randomly,  the  probability  of  detecting  a  user  when  k  tests  are  administered  to 
him  is  given  by 

This  is  the  probability  of  at  least  one  detection  among  k  tests.  To  minimize  the 
cost  of  testing  to  meet  the  requirement  of  identifying  lOOf  percent  of  the  users, 
the  smallest  k  for  which 

l-d-Dul* 

is  selected.  Then,  k  is  the  smallest  nonnegative  integer  satisf3ring 


Note  that  k  can  be  viewed  as  the  number  of  tests  an  individual  is  to  be  given 
over  his  lifetime  in  the  service.  In  actuality,  these  tests  would  be  given  within 
a  specified  period  (for  example,  a  year),  during  which  they  could  be  admin¬ 
istered  in  a  random  fashion  to  a  relatively  stable  population. 


Example  1 

The  following  example  illustrates  the  use  of  the  model  in  establishing  a 
testing  program  at  minimal  cost.  For  this  illustration,  aribitrary  values  are 
assigned  to  the  variables. 

The  testing  population  is  divided  into  five  classes.  The  distribution  of 
usage  classes  for  the  population  is  assumed  to  be  known  and  to  be  as  follows; 


Class  Usage  rate  Percent  of  population 

0  0  (nonusers)  75 

1  Once  per  month  10 

2  Once  per  week  8 

3  Twice  per  week  5 

4  Once  every  2  days  2 


This  information  can  be  used  as  follows  to  estimate  the  parameters  of  the 


model; 

Class 

0) 

Probability  of  use 
on  arbitrary  day 

pO) 

Percent  of 
population 

fiJ) 

Percent  of  user 
population 
SO) 

0 

0 

75 

0 

1 

12/356=  0.03 

10 

40 

2 

1/7  =  0.14 

8 

32 

3 

2/7  =  0.29 

5 

20 

4 

1/2  =  0.50 

2 

8 

Note  that  we  have  interpreted  the  usage  rate  as  the  probability  that  a  drug 
will  be  used  on  a  randomly  chosen  day. 


For  this  example,  the  test  is  assumed  to  be  effective  with  a  probability 
of  1  if  the  person  being  tested  last  used  the  drug  3  days  at  the  most  before  the 
test  is  given  to  him. 

Table  1  gives  the  steady-state  probabilities,  x(i),  and  the  detection  rate 
for  each  class. 


TABLE  1 

STEADY-STATE  PROBABIUT1ES  FOR  STATES 
AND  DETECTION  RATES 


States 


Class 

1 

x(1) 

2 

x(2) 

3 

x(3) 

4 

x(4) 

Probability 
of  detection 
DO) 

0 

0.00 

■1 

1.00 

0.00 

1 

0.03 

mgsm 

0.03 

0.91 

0.09 

2 

0.14 

0.12 

0.10 

0.64 

0.36 

3 

0.29 

0.20 

0.15 

0.36 

0.64 

4 

0.50 

0.25 

0.125 

0.125 

0.88 

The  detection  rates  for  the  entire  population  (D)  and  the  user  population 
(Du)  for  a  single  test  are 

D  =  0.09 
Du  =  0.35  . 

The  desired  number  of  tests  per  individual  to  identify  95  percent  of  users  is  the 
smallest  nonnegative  integer  k  for  which 


.  /n(l-0.95) 

k  a  ■: -  =  6.95 

Wl-0.35) 


Hence,  we  take  fe  =  7. 
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Example  2 


In  the  previous  example,  it  was  assumed  that  usage  on  any  given  day  is 
described  by  a  Bernoulli  distribution  whose  mean  for  class  j  is  the  usage  rate. 
It  is  interesting  to  compare  the  above  results  with  those  obtained  when  drug 
usage  is  assximed  to  follow  another  distribution.  If  it  is  assumed  that  usage 
occurs  randomly  (uniformly)  over  time,  the  time  between  consecutive  uses  is 
exponential.  The  mean  of  the  exponential  distribution  describing  the  usage 
for  class  j,j  >  0,  will  be  the  reciprocal  of  the  usage  rate  for  class  j  .  Under  these 
assumptions,  the  probabilities  of  use  on  an  arbitrary  day  are  the  following; 


Class 

0) 


Probability  of  use 
on  arbitrary  day 

pO) 


Percent  of 
population 

r(i) 


Percent  of  user 
population 

s(j) 


0 

1 

2 

3 

4 


0.00 

0.03 

0.13 

0.25 

0.39 


75 

10 

8 

5 

2 


0 

40 

32 

20 

8 


The  steady-state  probabilities  and  the  detection  rates  are  given  in 
table  2. 


TABLE  2 

STEADY-STATE  PROBABILITIES  FOR  STATES  AND  DETECTION  RATES 

States 


Class 

1 

x(1) 

2 

x(2) 

3 

x(3) 

4 

x(4) 

Probability  of 
detection 

OO) 

0 

0.00 

0.00 

0.00 

1.00 

0.00 

1 

0.03 

0.03 

0.03 

0.91 

0.09 

2 

0.14 

0.12 

0.10 

0.64 

0.36 

3 

0.25 

0.19 

0.14 

0.42 

0.58 

4 

0.39 

0.24 

0.15 

0.22 

0.78 
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The  detection  rates  for  the  entire  population  (D)  and  the  user  population  (Du) 


D  =  0.08 


Du  =  0.34  . 


The  minimum  number  of  tests  per  individual  needed  to  detect  95  percent 
of  the  actual  users  is  8. 

APPLICATIONS  OF  THE  MODEL 

Three  additional  applications  of  the  Markov  chain  model  are  presented 
below.  The  applications  are  directed  at  determining  the  percentage  of  drug 
users  in  the  population,  estimating  the  particular  usage  class  to  which  an 
indivdiual  belongs,  and  estimating  the  steady-state  population  in  each  usage 
class  when  a  policy  of  discharging  detected  users  is  followed. 

Estimating  the  Percentage  of  Drug  Users  in  a  Population 

Drug  testing  is  a  deterrent  to  drug  usage.  Historical  evidence  suggests 
that  when  testing  is  used  and  detected  drug  users  are  penalized,  the  decrease 
in  the  percentage  of  the  population  in  a  particular  class  of  users,  other  than 
the  class  of  nonusers,  is  independent  of  the  class.  That  is,  if  the  population  of 
class  j,j  >  0,  decreases  by  x  percent,  then  the  population  in  every  class  other 
than  0  will  also  decrease  by  x  percent. 

Assiune  that  the  rij)  values  represent  the  decimal  equivalent  of  the 
percentage  of  the  population  in  each  class  before  the  testing  program  begins. 
After  the  program  has  begun,  the  percentages  in  each  class  are  expected  to 
change.  Let  t(j),j  =  0,  1,  ...,  m,  denote  the  portion  of  the  population  in  class  j 
after  program  initiation.  The  assumption  of  proportionality  implies  that 
there  is  a  constant  c  for  which 


Hjl  =  ci<j)  0  =  1.2 . ffi) 


«0)  =  l-«(l)  +  «(2)  +  ... +  <(m)l  . 


If  D  is  the  detection  rate  when  the  distribution  of  users  is  given  by  rij) 
and  Do  is  the  detection  rate  when  the  distribution  of  users  is  given  by  t(j),  then 

Do  =  <(0)D(0)  +  /(l)D(l)  +  ...  +  Km)D(m) 

=  dmOiO)  +  r(l)7>(l)  +  ...  +  Km)  Dim)] 

=  cD  , 


Example  3  (Example  I  continued) 

Suppose  that  after  some  time  has  passed  since  testing  began,  the 
detection  rate  (Do)  is  found  to  be  0.06  instead  of  the  predicted  (D)  0.08.  Then, 


0.06 

c=  - =  0.75  . 

0.08 

Using  this  value,  it  is  possible  to  estimate  the  percentage  of  the  population  in 
each  usage  class.  Table  3  gives  the  original  percentages  and  the  percentages 
for  various  detection  rates  arrived  at  over  time  after  testing  began. 


TABLE  3 

PERCENTAGE  OF  POPULATION  IN  CLASS  /  FOR  VARIOUS  DETECTION  RATES 

Detection  rate 


p(i) 

0.08 

0.06 

0.04 

0.02 

0 

75 

81.2 

88.0 

93.4 

1/2 

10 

7.5 

5.0 

2.4 

1/7 

8 

6.0 

4.0 

2.0 

‘.•V, 

2/7 

5 

3.8 

2.5 

1.3 

12/365 

2 

1.5 

1.0 

0.5 

^  A 


Estimating  the  Usage  Class  of  an  Individual  Based  on  Testing 


Assume  that  k  tests  will  be  given  to  each  individual  in  the  population. 
How  can  the  usage  class  to  which  an  individual  belongs  be  estimated  based  on 
the  outcome  of  these  tests?  If  an  individual  has  been  given  g  rs  k  tests  and 
Q^f^gof  them  have  been  positive,  what  can  be  concluded  about  the  class  to 
which  this  individual  belongs? 

The  probability  of  an  individual  being  in  class  given  that  font  of  g  tests 
given  to  him  have  been  positive,  P(/),  can  be  determined  as  follows: 


P(/)  =  P((n  class  j  given  fofg  tests  positive) 

P(ifi  class  j  and  fofg  tests  positive) 
Pifofg  tests  positive) 

_  r</)C(g.fl[D(//m-DO)f^-^ 
a 

“  ~  * 

w 


where 


u  =  r(0)C(g,flD(0/[l-Dm^-^  +  ...  +  Hn)C(g,flV(n/  [1-D(n)l^-^ 

and 

w  =  riO)  [1  -D(0)^-P  +  ...  +  iin)D(n/ [1  -0(n)]^-^  j  , 
where  C(g,  f)  denotes  the  binomial  coefficient  "g  binomial  f.” 

Example  4 

Assume  the  parameters  defined  in  example  1,  which  are  repeated  in  the 
following  table. 


Class 

01 

Probability  of  use 
on  arbitrary  day 

poi 

Percent  of 
population 

d/) 

Percent  of  user 
population 

SOI 

0 

0 

75 

0 

1 

12/365  =  0.03 

10 

40 

2 

1/7  =  0.14 

8 

32 

3 

2/7  =  0.29 

5 

20 

4 

1/2  =0.50 

2 

8 

The  detection  rates  are  as  follows; 

States 

1  2  3 

Class  x(1)  x(2)  je(3) 

4 

d4) 

Probability  of 
detection 

oo) 

0 

0.00 

0.00 

0.00 

1.00 

0.00 

1 

0.03 

0.03 

0.03 

0,91 

0.09 

2 

0.14 

0.12 

0.10 

0.64 

0.36 

3 

0.29 

0.20 

0.15 

0.36 

0.64 

4 

0.50 

0.25 

0.125 

0,125 

088 

Table  4  gives,  for  each  class,  the  probabilities  of  an  individual  being  in 
that  class,  based  on  1  to  4  tests  and  the  specific  number  of  positive  tests. 

Policy  Decisions 

In  using  the  model  to  estimate  the  percentage  of  drug  users  in  a 
population,  it  was  assumed  that  the  punishment  policy  used  to  reprimand 
detected  drug  users  resulted  in  a  distribution  of  users  among  classes  that 
maintained  a  fixed  proportionality  among  the  numbers  in  the  various  usage 
classes.  The  model  can  also  be  used  to  investigate  how  alternative  policy 
decisions,  where  proportionality  will  not  be  maintained,  will  affect  the 
distribution  of  the  entire  population  across  classes.  In  particular,  the  effect  of 
discharging  users  who  have  had  a  specified  number  of  positive  tests  within  a 
given  testing  period  is  examined. 


TABLE  4 


PROBABIUTY  OF  AN  INDIVIDUAL  BEING  IN  CLASS  0  TO  4,  BASED  ON  TESTING 


Probability 

Number  of  Number  of 


tests 

positive  tests 

Class  0 

Class  1 

Class  2 

Class  3 

Class  4 

1 

0 

0822 

0.100 

0.056 

0020 

0.003 

1 

0.000 

0.103 

0.330 

0.336 

0.201 

2 

0 

0860 

0.095 

0.038 

0.007 

0.000 

1 

0.000 

0.203 

0.458 

0.286 

0.052 

2 

0.000 

0.017 

0.220 

0434 

0.329 

3 

0 

0.884 

0.089 

0.025 

0.003 

0.000 

1 

0.000 

0.315 

0.499 

0.175 

0.011 

2 

0.000 

0.044 

0.400 

0.444 

0.112 

3 

0.000 

0.002 

0.122 

0.429 

0.446 

4 

0 

0.901 

0.082 

0.016 

0.001 

0.000 

1 

0.000 

0.428 

0.476 

0.094 

0.002 

2 

0.000 

0.086 

0.545 

0.341 

0.029 

3 

0.000 

0.008 

0.271 

0.536 

0.186 

4 

0.000 

0.000 

0.062 

0.386 

0.552 

Assume  that  the  distribution  of  drug  users  across  the  classes  of  people 
entering  the  Navy  from  the  civilian  population  is  known.  Let  Hj),J  =  0, m 
denote  the  probability  that  a  new  inductee  will  be  in  user  class  j. 

Let  w  represent  the  probability  that  an  arbitrarily  selected  individual 
will  leave  the  service  during  the  period  being  used  for  reasons  other  than  the 
drug  detection  policy  described  below. 

Also  asstune  that  in  any  given  period,  k  tests  will  be  administered  and 
that  any  person  detected  h  out  of  k,  h  ^  k,  times  in  this  period  will  be 
discharged  from  the  Navy. 


.12- 


The  probability  of  detection  within  each  class  can  be  used  to  compute  the 
probability  of  having  at  least  h  out  of  k  tests  positive  for  a  randomly  selected 
individual  in  that  class.  This  probability  is  denoted  by  P(j),  j  =  0, m. 

If  rij)J  =  0, m,  denotes  the  percentages  of  the  population  in  classes 
0, m,  respectively  at  steady  state  (i.e.,  after  the  policy  has  been  put  into 
effect  and  the  classes  have  stablized),  then  the  rate  at  which  individuals  will 
leave  class  j  is  given  by 


m  •  ■  •  •  • 


•  -  O 


•■.--•.'•V 


‘.•V-.i 


P{j)  rij)  +  «i;[l  -P(j)]  rij)  =  {PO")  +  u»[l  -P(;')l)d/')  =  B(j)  rij) 


where 


B{f)  =  P<j)+w[l-P{j)]  . 


The  rate  at  which  individuals  will  leave  the  service  is  given  by 


X  =  B(O)KO)  +  Ba)Kl)  +  ...  +  B(m)dm)  . 


The  rate  at  which  people  leave  the  service  is  assumed  to  be  equal  to  the 
rate  at  which  they  enter,  so  that  the  service  population  is  constant.  Since  it  is 
assumed  that  rij)  represents  the  steady-state  percentages  in  each  usage  class, 
the  rate  of  entry  into  each  usage  class  must  be  equal  to  the  exit  rate  for  each 
usage  class.  That  is. 


*  ■  •  W 

•  •  I 


tO')X  =  BO‘)K/l 


This  set  of  m  equations  for  the  m  -J-  1  unknowns  r{j),j  =  0, ...,  m  has  one 
redundant  equation.  However,  rij)  must  also  satisfy  the  requirement  that 

KO)  +  Kl)  +  ...  +  Km)  =  I  , 


V  V/." 

•.*  'V 

•  •  *  •  -  ' 
•  -J 

N 

‘"s-'VN 

fv-v* 


which  leads  to  the  following  system  of  equations: 


W0)-1]B(0)K0)  +  <(0)B(l)r<l)  +  ...  +  <(0)B(m-l)Km-l)  +  Hl)B{m)r{m)  =  0 
<(l)B(0)d0)  +  [<(l)-llB(l)f<l)  +  ...  +  <(l)B(m-l)Km-l)  +  <(l)B(m)Km)  =  0 


<(/n-l)B(0)/<0)  +  <(m-l)B(l)/<l)  +  ...  +  Wm-l)-llB(ni-l)/<m-l) 
+  Kn—  l)B(m)/<m)  =  0 

KO)  +  Kl)  +  Km— 1)  +  Km)  =  1  . 


Let 


A0‘)  =  B(O)...B(;-l)B0+l)...B(m)  0  =  1 . m) 


Then  the  solution  of  this  system  is  given  by 


K/)  = 


<01  AO) 


[<(0)A(0)  +  ...  +  <(m)A(m)] 


0  =  0,  ...,m) 


Example  5  (Example  1  continued) 

The  rate  at  which  people  leave  the  population  for  reasons  not  related  to 
drug  detection  is  assumed  to  be  u;  =  0.1  (i.e.,  10  percent  of  the  population  will 
leave  if  no  new  drug  policy  is  followed).  It  is  also  assumed  that  new  inductees 


obey  the  distribution  of  drug  users  described  in  example  1,  which  is  as  follows: 


Probability  of  use 
on  arbitrary  day 

pO) 


12/365=  0.03 
1/7  =  0.14 

n 

12 


Percent  of  inductee 
population 


Percent  of  inductee 
user  population 


The  following  table  gives  the  detection  rates  for  the  usage  classes  for  a 
single  test: 


States 


Class 

1 

x(1) 

2 

x(2) 

3 

x(3) 

4 

x(4) 

0 

0.00 

0.00 

0.00 

1.00 

1 

0.03 

0.03 

0.03 

0.91 

2 

0.14 

0.12 

0.10 

0.64 

3 

0.29 

0.20 

0.15 

0.36 

4 

0.50 

0.25 

0.125 

0.125 

Probability  of 
detection 

m 


••  .•  V 


Now  suppose  that  seven  tests  are  given  at  random  over  a  period  of  one 
year,  and  the  policy  is  to  discharge  anyone  who  is  detected  on  at  least  two  of 
the  tests.  The  probabilities  ofdischarge  for  each  class,  P0)>  are 


P(0)  =  0.000 


P(l)  =  0.125 


P(2)  =  0.783 


P(3)  =  0.826 


P(4)  =  l.OOO  . 


».•  V*  •  J 


."k , 


Note  that 


fl(0)  =  0  +  (l-OKO.l)  =0.1000 

B(l)  =  0.125  +  (1  -0.125)(0.1)  =  0.2125 
B(2)  =  0.8047 

B(3)  =  0.8434 

B(4)  =  0.1000 


A(0) 

A(l) 

A(2) 

A{3) 

A(4) 


(0.2125)(0.8047X0.8434)(1)  =  0.1442 
(0.1000)(0.8047K0.8434)(1)  =  0.0679 

=  0.0179 
=  0.0171 
7  0.0144 


and 


4(0)A(0)  =  (0.75X0.1442)  =  0.1082 
<(1)A(1)  =  (0.10X0.0679)  =  0.0068 
<(2)A(2)  =0.0014 

t(3)A(3)  =0.0009 

<(4)A(4)  =0.0003 

/(O)A(O)  +  ...  +  <(4)A(4)  =  0.1175  . 


Hence,  the  steady-state  percentages  in  the  usage  classes  are 


liO)  =  0.1082/0.1175  =  0.921 
Kl)  =  0.0068/0.1175  =  0.058 
K2)  =0.012 

r<3)  =  0.007 

K4)  =  0.003  . 


Summary 

These  applications  of  the  Markov  model  illustrate  the  potential  of  the 
model  for  analyzing  a  population  in  terms  of  drug  usage.  Given  a  drug- testing 
program,  the  model  can  be  used  to  estimate  the  distribution  of  drug  users  in 
the  population  on  the  basis  of  observed  detection  rates  and  to  identify  the 
usage  class  of  individuals  based  on  testing  results.  Finally,  the  last  example 
shows  how  the  model  can  be  used  to  predict  the  effects  of  policy  decisions  on 
drug  use. 
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